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Motivated by a long-standing debate concerning the nature and interrelations of surface-tension 
variables in fluid membranes, we reformulate the thermodynamics of a membrane vesicle as a generic 
two-dimensional finite system enclosing a three-dimensional volume. The formulation is shown 
to require two tension variables, conjugate to the intensive constraints of area per molecule and 
volume-to-area ratio. We obtain the relation between these two variables in various scenarios, 
as well as their correspondence to other definitions of tension variables for membranes. Several 
controversies related to membrane tension are thereby resolved on a model-free thermodynamic 
level. The thermodynamic formulation may be useful also for treating large-scale properties of 
vesicles that are insensitive to the membrane's detailed statistical mechanics and interactions. 
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Membrane vesicles are flexible fluid envelopes, made of a bilayer of amphiphilic molecules in aqueous solution [H, 0. 
Since Helfrich's work almost three decades ago Q , these prevalent structures have been one of the most extensively 
studied systems in soft-matter physics, exhibiting rich statistical-mechanical [|J and dynamic Q behaviors. It may 
seem odd, therefore, to go back now and reformulate their equilibrium thermodynamics. Nevertheless, this is the goal 
of the current work. The reason is a long and ongoing debate concerning various definitions of the surface tension(s) of 
fluid membranes, their physical meanings, and interrelations [6H20j. Such issues as the exact number of independent 
tension variables and their possible vanishing ((J Qjl, [HI , or even sign reversal > have remained controversial. 

The debate has been conducted at the statistical-mechanical level, emphasizing the role of thermal fluctuations and 
employing various analytical and simulation techniques. While certain issues are inherently statistical-mechanical, it 
seems plausible that at least some of the other fundamental issues mentioned above should be settled at the lower, 
more general level of thermodynamics (i.e., without specifying a Hamiltonian) . 

Two key properties of fluid membranes make their thermodynamic characterization special, and confusing. The first 
is that, while being extensive in two dimensions, the membrane is embedded in three dimensions. Consequently, it 
was suggested early on that the thermodynamics of open membranes required two spatial constraints rather than one 
| 6j — the membrane area, A, and its projection onto a reference planar frame, A p — entailing two conjugate variables 
8J. The variable conjugate to A, II, is a surface pressure, analogous to ordinary thermodynamic pressure. The one 
£S) " related to A p , the frame tension r, is the force per unit length exerted at the edges of the frame, in the plane of 
reference, to maintain a given A p . 

The second property of bilayer membranes which sets them apart from ordinary fluids is the fact that they have a 
y—t \ preferred area per molecule [l2|, |21| . This boils down to the hydrophobic effect, which opposes the exposure of the 



y—i ■ inner core of the bilayer to the outer solvent [21|. As a result, the membrane possesses in-plane elasticity — unlike an 
ordinary fluid it can sustain not only a global compressive stress but also a global dilative ones. 

• i-H ■ The discussion of membrane thermodynamics is further obscured when the in-plane elasticity is suppressed by 
" taking the limit of membrane incompressibility. This commonly used assumption simplifies calculations and allows 
one to theoretically fix the total area of the membrane. Yet, at the same time, it makes essentially independent 
thermodynamic variables dependent — the area A and number of membrane molecules N become linearly related, 
and so do the conjugate surface pressure and chemical potential (3, Il2j. 

To these complications, which are special to fluid membranes, one should add the common pitfalls in thermodynamic 
analyses — the necessity to precisely define which variables are kept fixed when another variable is varied, and the 
fact that two variables, which turn out to be equal under certain considerations, are not necessarily identical. 

To circumvent some of the foregoing complications we consider a closed envelope — a vesicle — rather than an open 
membrane. We expect (and will demonstrate below) that, once the thermodynamic behavior of a vesicle is clarified, 
the consequences for an open membrane should be readily inferred. In the case of a closed vesicle the three-dimensional 
constraint is an ordinary volume constraint on the enclosed system, rather than a projected-area constraint on the 
surface. Thus, we explicitly consider the contents of the envelope as a thermodynamic system, for which the properties 
of the membrane become surface effects. As a result, the starting point of the formulation resembles that of other 
studies of finite-size systems such as droplets or domains. It contains only well known, unambiguous thermodynamic 
variables, such as volume and surface area. The special properties of the membrane subsequently emerge without 
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further assumptions or an externally imposed reference frame. In addition, membrane incompressibility is nowhere 
assumed. 

The current work should be distinguished from other thermodynamic formulations for membranes, which were 
presented in the past (22l . [23l. The aim of those theories was to use a local thermodynamic description, based on 
Gibbs' theory of interfaces [22[ or molecular considerations 23 1, to construct a coarse-grained, local free energy of a 
non-fluctuating mesoscopic element of the membrane. Subsequently, the resulting free energy per unit area could be 
used, like the Helfrich Hamiltonian Q , as a starting point for statistical- mechanical or other calculations on a larger 
scale. By contrast, we directly address the global, large-scale thermodynamics of a membrane vesicle. 

In Sec. [IT] we formulate the thermodynamics of the vesicle and present the resulting tension variables. In Sec. Mil 
we examine the consequences of the formulation in various experimental and theoretical scenarios, and the resulting 
relations between the two tension variables themselves, and between them and other previously defined variables. 
While the formulation is for a closed vesicle, we examine the analogous results for open membranes as well. Section 
IIVI summarizes the findings, relates them with earlier results, and discusses their implications. 



II. THERMODYNAMIC FORMULATION 

Consider a three-dimensional system of volume V , containing Q particles 24] and being in thermal contact with a 
bath of temperature T. In the thermodynamic limit the Helmholtz free energy of the system, F 3 (T, V, Q), is extensive 
in Q while keeping the volume per particle fixed, 

F 3 (T,V,Q)=Qf 3 (T,v), v = V/Q. (1) 

If we double the number of particles, and at the same time also double the volume, the free energy F% will double. 
Conjugate to T is the entropy per particle, S3 = —{df 3 /dT) v . Conjugate to v is the internal pressure, 

The volume is enclosed by a two-dimensional envelope of surface area A and N particles [25[ . As in studies of other 
finite-size systems (e.g., droplets, bubbles, domains), we restrict the analysis to the leading correction introduced 
by the surface to the thermodynamic limit F — F3 of Eq. ([l}. This leading correction arises from contributions 
to the free energy which are extensive in the surface size N. Since the envelope is closed in three dimensions, its 
Helmholtz free energy generally depends on both A and V, F 2 — F 2 (T, A, V, N). For example, the envelope's entropy 
(which is extensive in N and should be included, therefore, in F 2 ) will decrease if V is increased while keeping A 
fixed. The dependence of F 2 on A and V implies that the thermodynamic description of the envelope requires two 
intensive variables apart from temperature. We choose these variables to be the area per particle, a = A/N, and 
the dimensionless volume-to-area ratio, a = 6^/ttV/A 3 ^ 2 (defined such that it has the maximum value of unity for a 
perfectly spherical envelope). For F 2 to be extensive in N, one should fix both a and a, 

F 2 (T,A,V,N)=Nf 2 (T,a,a), a = A/N, a = 6^V/A 3 / 2 . (3) 

This is a key point which is worth dwelling upon. Imagine, for example, that we double both the number of particles 
making the envelope and its area, but do not change the enclosed volume. The surface free energy arising from 
in-plane short-range interactions and in-plane entropy will double. Yet, evidently, the resulting object is not merely 
a rescaled version of the original one. It is less spherical and bound to have more out-of-plane fluctuation entropy 
per molecule, leading to more than double the surface free energy. If we double N and A, and at the same time 
also increase V by a factor of 2 3 / 2 , the envelope will be properly rescaled. The in-plane as well as the out-of-planc 
contributions per molecule will remain unchanged, and the surface free energy F 2 will double. We stress that Eq. Q 
contains only contributions which are extensive in the envelope's size. 

The appearance of the thermodynamic variable a reflects a three-dimensional constraint imposed on the envelope. 
Here the constraint is simply that the envelope be closed and contain a volume V. One may impose other three- 
dimensional constraints, such as a given projected area, A p , onto a reference surface. Such a thermodynamic constraint 
is not independent of the volume constraint; it may replace it but not be added to it. For example, if the area A of 
the envelope and the volume V that it encloses are known, then, evidently, the area of the sphere of volume V (taken 
as a reference surface), A p = (47r) 1 / 3 (3^) 2 ^ 3 , and its deviation from A, are known as well. It is clear, however, that 
the foregoing considerations of A and V are more transparent and independent of any reference frame. 

The flexibility of the thermodynamic formulation allows for adding other constraints. For example, one may suggest 
that F 2 should depend on additional characteristics of membrane shape apart from A and V, such as the (average) 
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curvature. Yet, one should examine whether such additional constraints are beneficial or experimentally relevant. It 
is hard to conceive a practical scenario where the shape of a membrane is controlled to such an extent. We prefer, 
therefore, to remain with the standard constraints of volume and area as considered for any finite-size system. In 
statistical-mechanical terms the membrane will sample all those additional shape characteristics which are consistent 
with the global constraints of A and V (with appropriate Boltzmann weights). 

Returning to Eq. ([3]), we write the entropy per particle of the envelope as s 2 = —(df2/dT) a ,of The dependence 
of f 2 on two intensive variables other than temperature, a and a, implies that there are necessarily two independent 
conjugate variables: 



II(T, a, a) — — 
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(4) 
(5) 



The prefactor in the definition of 7 has been introduced to make the following presentation clearer. Both variables have 
surface-tension dimensions. The first, II, is the surface pressure — the two-dimensional analogue of p m — describing 
the response of the envelope's molecules to in-plane strain (i.e., changes in the area per molecule without a change in 
the volume-to-area ratio). At the optimum area per molecule II vanishes. The second variable, 7, referred to hereafter 
as the Laplace tension, is related to the response of the envelope to out-of- plane strain (i.e., changes in the enclosed 
volume without a change in the area per molecule). 
We can further define three surface moduli, 
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of which only the first is usually considered (aK aa is the membrane's compression modulus). The last equality 
in Eq. ([6]) is a Maxwell relation between the two variables, II and 7. In addition, convexity of f 2 demands that 
K K > K 2 

Now consider the combined system, with total Helmholtz free energy F(T, V, Q, A, N). The fundamental differential 
relation for F is 



dF = -SdT - pdV + n 3 dQ - odA + fi 2 dN, 



(7) 



where S is the total entropy, p the total pressure exerted by the environment, a the total surface pressure in the 
envelope, and fj,^ and /12 the chemical potentials, respectively, of the particles enclosed inside the vesicle and those 
making the envelope. In what follows we assume for simplicity that the particles in the volume interact with those 
in the envelope via hard-core repulsion only, i.e., they are merely enclosed by the envelope. In this case the total 
Helmholtz free energy is simply given by the sum of the volume and surface contributions, 



F = F 3 +F 2 = Qf 3 (T, v) + Nf 2 (T, a, a). 



(8) 



Note, however, that there is an implicit coupling between the two terms via their mutual dependence on the volume. 
Taking the differential of Eq. ([5]) and equating it with Eq. ([7]), we identify: 
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= Qs 3 + Ns 2 , 


(9) 




2 [NaY 1/2 


(10) 


P 


= P --a{-^) 7 ' 




= f3+PinV, 


(11) 


a 


- n + 7, 


(12) 


/'2 


= f2 + Ua. 


(13) 



Thus, the total surface pressure a — the variable conjugate to A in the total free energy — is manifestly broken 
into two distinct contributions, II and 7 [Eq. (IT21) ]. We underline the different roles played by each of these two 
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tension-like variables, denned in Eqs. ((4]) and (0). The similarity between II and ordinary thermodynamic pressure 
is reflected in its relation to the chemical potential of the envelope's molecules, Eq. (fT5|). in which 7 is absent. The 
Gibbs free energy of the envelope alone, whose sole extensive variable is N, is G2 = H2N = F2 + HA (not F2 + a A). 
By contrast, the mechanical nature of 7 is reflected in its contribution to the pressure that the vesicle exerts on the 
external environment, Eq. (|10p . in which II does not appear. We may define for the vesicle an effective thermodynamic 
radius of curvature, as 

fNa\ 1/2 W 

which for a perfectly spherical envelope coincides with the sphere's radius. With this definition Eq. (|TU|) becomes 
Laplace's law, p m — p = 2j/R c . 

III. EXPERIMENTAL AND THEORETICAL SCENARIOS 

As in any thermodynamic analysis, we can apply the formulation given above to various scenarios, in which certain 
variables are controlled while others are free to relax. These scenarios correspond to different statistical ensembles. 
In the following subsections we describe several scenarios of particular relevance to experimental conditions or to the 
debate over membrane tension. 



A. Volume or pressure constraint 

Let us begin with a scenario in which we prescribe the volume enclosed by the envelope, as well as the temperature 
and the number of particles in the volume and in the envelope — i.e., we fix (T, V, Q, N). In practice the volume may 
be considered fixed if we assume the inner solvent to be incompressible and neglect the permeation of solvent into and 
out of the vesicle. (This is valid, for example, over sufficiently short time.) Previous studies considered an additional 
constraint of either a fixed area A or a fixed total surface pressure a. We argue that this is unnecessary. Indeed, it is 
impractical to dictate the area of the membrane. If we do not impose an external surface pressure (a = 0; the case of 
finite a is discussed in the next subsection), the area will relax under the given constraints of (T, V, Q, N) 26]. The 
Helmholtz free energy is then minimized with respect to A, leading to 

n + 7 = o. (15) 

After the minimization, F remains dependent only on (T, V, Q, N). The force balance expressed by Eq. (|T5|) has the 
following intuitive interpretation. The volume constraint acts as if an external surface pressure were applied to the 
envelope, n ox = —7 (27j . In response, the particles of the envelope develop an internal surface pressure, II = n cx . The 
membrane's response to out-of-plane strains must be restoring, 7 > 0. Hence, in the absence of an actual external 
surface pressure, the in-plane balance is between two negative pressures — i.e., the volume constraint acts to increase 
the area, whereas the internal pressure acts to decrease it. 

Note that Eq. (|15[) is not an identity but an equation of state, relating the envelope's intensive variables, T, a, and 
a: II and 7 are in general independent variables, which become related in this scenario due to area relaxation. In 
addition, using Eqs. (|5J), @, jSj, and ([15]). together with the known dependence of a on V and A, one can calculate 
for a given model the unknown variables, such as p, a, II, and 7, as functions of the constraints (T, V, Q, N) [28| . 

Now consider a different case, where the pressure of the external environment is controlled, while the volume can 
vary, i.e., we fix (T,p, Q, N). This scenario is quite similar to the preceding one; we discuss it separately because it is 
the most relevant for actual equilibrium vesicles, where the solvent has sufficient time to permeate in and out of the 
vesicle. In this case we should minimize the Gibbs free energy, G = F + pV, with respect to V and A, leading again 
to Eqs. I|10p and l|15p. After the minimization G remains dependent only on (T,p, Q, N). We note that, although we 
dictate p, we do not directly fix p ln because of the surface effect represented by the Laplace tension 7 in Eq. (fTU)) . 

If there is neither a volume constraint nor a pressure one (i.e., p — Pi n ), we have a tensionless membrane, II = 7 = 

1- 

B. Contact with a reservoir of amphiphilic molecules 

Let us examine a scenario in which the set of constraints is (T, p, Q, ^2)- We explicitly consider this case for two 
reasons. First, it becomes the practically relevant scenario at sufficiently long times, when the vesicle exchanges 
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amphiphilic molecules with the surrounding solution. Second, it highlights the crucial role played by the three- 
dimensional embedding of the membrane. The latter issue is recognized once we notice that the set (T,p,Q, n 2 ) does 
not include any constraint that is extensive in membrane size. On the one hand, such a thermodynamic formulation 
would normally be under-defined and, hence, invalid. In the absence of any constraint on the envelope's size, the 
vesicle can grow indefinitely by drawing more amphiphilic molecules from the reservoir. On the other hand, it is clear 
that this scenario can be realized in experiments and simulations. This apparent contradiction is resolved below. 

The free energy to minimize in this case is G = F + pV — ^iV . Minimization of G with respect to V, A, and N 
recovers Eqs. (fTU|) . p^|) . and (IT51) . which we should be able to solve for V, A, and N, given the constraints (T,p, Q, fj,^)- 
Let us explicitly rewrite the equations in terms of these unknowns: 

° = h { T 4^) + ^ n { T 4^) - ^ 
° = n { t 4>Jji) ^ { T 4>Jn) > 

f V\ 2A ( A V \ 

The first two equations, describing in-plane equilibrium, are invariant to the scaling (N — > XN, A — > XA, V — > X 3 ^ 2 V) 
by an arbitrary scale factor A. It is the third equation (Laplace's law), accounting for the out-of-plane force balance, 
which violates this scaling and prevents the set of equations from being under-defined. 



C. Surface pressure constraint 



Just as it is impractical to impose a fixed area A on a vesicle, it is unclear how one could in practice apply to 
the vesicle a fixed in-plane force per unit length, a [2j|. Nonetheless, in statistical-mechanical studies of membranes 
such an external surface-pressure term, a A, is frequently added to the Hamiltonian as a means to control the mean 
area. In the current discussion there is no need to introduce such a term and, moreover, we shall see that it leads 
to unnecessary complications; we consider it merely to relate the emergent variables, LI and 7, with the commonly 
used a. Within our thermodynamic formulation, such a modification adds a term a A to the thermodynamic internal 
energy and free energy 30]. Consequently, minimization of the free energy with respect to A now yields Eq. (|12l) . 
a = II + 7, as the equation of state instead of Eq. (fT5|) . One can again interpret this result as if the particles of the 
envelope responded through an internal surface pressure LI to an external surface pressure, which consists of a and 
the pressure —7 due to the volume constraint, LI = II CX = a — 7. Thus, just as prescribing p in Sec. IIII Al did not fix 
Pin because of the existence of 7, introducing the external surface pressure a does not directly fix the internal one, LI, 
for a similar reason. 



D. Projected area and frame tension 



For open membranes a constraint on the membrane's projected area, A p , onto a certain reference plane, is widely 
used. In the case of a vesicle the choice of the reference surface itself is nontrivial and may depend on the specific 
scenario (la ]. We emphasize again that such a projection constraint is unnecessary in the current analysis since, as 
noted in Sec. [ill constraints on A p and V are mutually dependent. We may introduce a projection constraint instead 
of the volume constraint (on V or p) to clarify the relations between the current formulation and previous ones. 
Possible choices of a reference are the mean surface about which the vesicle fluctuates [HI, [3l[ , the spherical shell that 
has the same area as the vesicle [10| . or the spherical shell that holds the same volume as the vesicle [l(J [H, 32]. We 
choose the last option as the reference surface, whose area is 

A p = a 2 / 3 A. (16) 

We may also define the appropriate intensive variable, a p = A p /N . Lhus, fixing the projected area corresponds to 
fixing neither a nor a separately, but rather the product c?l 3 a. 

Related to the projected area is the frame tension, t. We have not used the term 'conjugate' intentionally, because 
the variable conjugate to A p would be the partial derivative of the membrane's free energy with respect to A p while 
keeping T, TV, and A fixed, 

ris (m =(m .(%>) (*>.) = a -v H (in 



dA p J TAN \da p J T a \ <9a / ^ „ \ <9a r 



6 



However, the frame tension as measured or controlled in experiments and simulations is related to the reversible work 
required to change A p while keeping only T and N fixed [13| . as it is impractical to fix A. This corresponds to 
differentiating fa with respect to a p while keeping only T fixed; both a and a are allowed to vary while satisfying Eq. 
(UHJ) and the equation of state a — a(T, a). (We have been assuming here the absence of an external surface pressure; 
if a nonzero a is considered, a should be kept constant as well.) This procedure yields 
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The general expression for r, Eq. (|18j) . contains a parameter e, which requires knowledge of the equation of state 
a(T, a). This parameter vanishes for an incompressible membrane having fixed a. It also cancels out in Eq. (|18p if we 
use the equation of state (|15l) (i.e., if we consider a = 0). In these two cases, therefore, we get the simple relation, 



a = or incompressible membrane: 



r = ti = a 2/3 7 = (A/A p )~f > 7. 



(19) 



Another simple limit is when the vesicle approaches a perfect spherical shape (or the open membrane approaches a 
flat state), where a — > 1 and e —¥ 00. Then we have r = —II, i.e., the frame tension provides the entire pressure 
required to balance the internal surface pressure. In the absence of an external surface pressure, a = 0, this gives the 
expected convergence for a = 1 of the frame and normal-response (Laplace) tensions, r = 7. 

If a nonzero a is considered, however, Eqs. (|12p and (|T5|) yield a different relation between the frame tension and 
the other tension variables, 



a- 2 / 3 [ 7 - 



1 



(20) 



which requires model-dependent knowledge of e. The parameter e can also be represented in terms of the surface 
moduli defined in Eq. ([6]). Using the equation of state (fT5|) or (fT2|) to write d(H + 7) = and extracting from it the 
derivative (da/da)T, we obtain 



3 (3/2)a 2 K aa — aaK aa + cry 
2 a 2 K aa - (3/2)aaK aa + 07 



(21) 



To assess the typical values of e we make the following estimates: we neglect the coupling modulus K aa ; for a vesicle 
approaching a spherical shape K aa ~ k^T j (1 — a) 2 [33|, where k^T is the thermal energy; 7a does not exceed a few 
k-^T before the membrane ruptures; K aa a 2 is of order 10 2 k-gT. Thus, for vesicles which are not extremely swollen, 
e <C 1, and the difference between Eq. (TT9l) for a = and ([20]) for a 7^ is small. For strongly swollen vesicles 
(1 — a < 10~ 2 ), however, e becomes large. In this case we get r — > 7 — a. Such a finite difference between the frame 
tension and the Laplace tension as the envelope and its reference surface coincide, is hard to physically justify; it is 
clearly related to the artificial inclusion of an external a 7^ 0. 



IV. DISCUSSION 



The requirement that a membrane vesicle, as complicated as its statistical mechanics may be, must comply with 
the thermodynamics of finite-size systems, has yielded a surprising amount of information. This information is in line 
with certain earlier findings and disagrees with others. 

The thermodynamics of a fluid vesicle requires two distinct tension variables, defined here as the surface pressure 
II and the Laplace tension 7 [Eqs. §4$ and ((5}]. We did not assume their existence but obtained them as a result of 
the standard geometrical constraints on the area and volume of the vesicle. When the area per molecule is allowed to 
relax — the relevant scenario in actual systems and molecular dynamics simulations — the total surface stress is zero, 
a = If + 7 = 0, and the two tension variables become related. Unlike early suggestions @, this does not imply that 
the membrane is tensionless, but just reflects the balance between the surface stress arising from a three-dimensional 
constraint and an equal internal surface pressure. As already recognized in Ref. Q, a state of vanishing tension, in 
the sense that the membrane does not experience a normal restoring force due to tension, is where 7 = 0. Then, upon 
area relaxation, we also have LT = 0, i.e., the area per molecule relaxes to its optimum value. This is achieved in a 
vesicle, for instance, when its volume is allowed to relax in conditions of zero pressure difference, p m = p [Eq. (|10[) ]. 
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A similar state is achieved in an open membrane (unless a a 7^ is introduced, or a non-optimum area per molecule 
is imposed), when its projected area is allowed to relax under conditions of zero frame tension, r = [Eq. Q19p], 

Another tension coefficient that is frequently used in the literature is the fluctuation tension, or '^-coefficient', r. It 
is extracted from the membrane's normal fluctuation spectrum, (u(q)u(— q)) = k-gT /[rq 2 + 0(q 4 )}, where u(q) is the 
Fourier-transformed normal displacement of the membrane as a function of wavevector q. Evidently, this definition 
of tension is purely statistical-mechanical and cannot be directly examined here. It is clear, nonetheless, that r is 
the tension related to the membrane's response to normal strain. We thus suggest to identify it as r = 7. Indeed, 
micropipette-aspiration experiments have shown that the measured r is consistent with the imposed vesicle's Laplace 
tension (33. For open membranes it was analytically argued I9I [L?! Il3l . IT9I l20l l , an d numerically demonstrated 
[T3I I20I . |35| . that r = t, the frame tension, but other works disagree [l4 . Il6l . l36l l37l |. Our analysis confirms the 
former results. We have found that, when there is no external surface pressure, a — 0, or when the membrane's 
compressibility can be neglected, the Laplace tension 7 is proportional to, but slightly smaller than, the frame tension 
t, 7 = a 2 / 3 r = (A p /A)t. This is in line with the reports of r — r, since the difference (r — 7)u(q)u(— q) is only 
fourth-order in the normal displacement. It suffices to redefine the Fourier transform using the membrane's material 
coordinates rather than its projected ones to get a slightly modified spectrum with a g 2 -coefficicnt r' = (A p /A)r = 
(A p /A)t = 7. A recent numerical study suggests, in fact, that the fluctuation tension is equal to (A p /A)t rather 
than t [L9| . As the vesicle becomes spherical (or an open membrane gets planar) the frame tension coincides with 
the Laplace tension, as expected. 

The term a A, added to the energy in many theories and simulations, is the origin of a lot of confusion. Thermo- 
dynamically, prescribing the variable a corresponds to applying an external surface pressure [30j |. similar to the pV 
term introduced when a three-dimensional system is in contact with a pressure bath. We are used to the fact that, 
for three-dimensional systems, we must prescribe either the volume or the external pressure; if we do not, there will 
be nothing to counterbalance the system's internal pressure. One initially expects that for a membrane we would 
similarly have to prescribe either A or a. (The latter is more widely used since studying fluctuations at fixed A is 
theoretically difficult [ToL l32l|.) However, in an actual vesicle nothing fixes the area, and there is no in-plane barrier 
that can apply a fixed external surface pressure (29[. Unlike the three-dimensional case, such an external pressure is 
not necessary to balance the internal pressure LI of the two-dimensional fluid, because the volume constraint (i.e., the 
embedding of the membrane in three dimensions) — be it the actual volume, a pressure difference, a projected area, 
or a frame tension — provides such a balancing pressure, equal to —7. Even if we tune the volume constraint down 
to zero (e.g., by imposing p nl = p or r = 0), the membrane can support a state of zero internal surface pressure by 
attaining its optimum area per molecule, without the need for an external pressure. Thus, the area is unconstrained, 
leading to IT + 7 = [28]; equivalently, we may say that the external surface pressure vanishes, a — 0. In this 
restricted sense, the thermodynamics of area changes is similar to the thermodynamics of photons, whose number is 
unconstrained or, equivalently, whose chemical potential vanishes. 

Hence, the ensembles of fixed (T,V,Q,N) or (T,p,Q,N) with a = (Sec. IIII Al) are the ones that we wish to 
advocate as the most physically relevant and the least confusing. The corresponding ensembles for open membranes 
are (T, A p , N) or (T, r, N), again with a = 0. (At sufficiently long times the constraint on N should be replaced by 
one on /i2', see Sec. IIII BH The distinction between ensembles is not merely technical or semantic, mainly because 
after removing the constraint on A or a the ensembles mentioned above contain one less constraint than the ones 
usually studied. In addition, even ensembles having the same number of constraints may not be equivalent for the 
finite systems under consideration, which contain both volume and surface free-energy contributions |19l . |38| . 

If one insists on including a nonzero a in the theory or simulation, then area relaxation will lead to a = II + 7, 
implying that a 7^ 7 = r in general. The inequalities a 7^ r and it/t were observed in simulations [HI dEH. We 
point out that the numerically measured difference, a — 7 = a — r, provides an indirect measurement of the internal 
surface pressure n in such nonzero-cr computations. The inclusion of a 7^ entails two additional problems which, as 
we have been trying to show, are superfluous, (a) Even when the membrane is tensionless, 7 = 0, there is a residual 
surface stress, II = a, causing the area per molecule not to relax to its optimum value, (b) The frame tension r is not 
proportional to 7 = r [Eq. (|2"0)) ]. Although the deviation should be usually small, when the vesicle is strongly swollen 
(or an open membrane is nearly flat), we get t — > 7 — <r, which is physically questionable for an envelope arbitrarily 
close to its projection. 

We end by commenting on the possible signs of the various tension variables [i!14i"8l |. The Laplace tension, which 
characterizes the membrane's normal restoring force, must be non-negative for the membrane to be stable, 7 > 0. 
This is reflected also in our identification 7 = r, since a negative r will cause unstable large-wavelength fluctuations. 
If we let the area relax without an external surface pressure, a = — the conditions which, as we have argued, are 
the valid ones for actual vesicles — then the condition 7 > entails a non-positive (compressive) internal pressure 
[Eq. (|15p], II < 0, and a non-negative frame tension [Eq. (|19p], r > 0. However, if a a 7^ is imposed, n can be also 
positive (dilative) [Eq. ([12"])]. and r can be also negative [Eq. ([2"0))]. 

The simple thermodynamic framework, laid out above, may be practically useful in cases where the detailed 
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statistical mechanics and interactions of a membranal system are not crucial. For example, a similar formulation has 
recently been used to unravel a universality in the osmotic swelling of vesicles, arising from the mere competition of 
volume and surface effects [!(| ■ There are issues which obviously cannot be dealt with using thermodynamics alone. 



and the equivalence (or lack of it) of different statistical ensembles [19|, [38|. Nonetheless, within its limitations, the 
remarkable strength of classical thermodynamics lies in its simplicity, broad applicability, and rigor [4l[. We hope 
that the formulation and conclusions presented here will be useful as guidelines for future studies. 
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discussions and correspondence with J.-B. Fournier, E. Haleva, M. Kozlov, J. Landy, and F. Schmid are gratefully 
acknowledged. Acknowledgment is made to the Donors of the American Chemical Society Petroleum Research Fund 
for partial support of this research (Grant No. 46748- AC6). 
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